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A new experimental evaluation of the quenching factor for nuclear recoils in NaI[Tl] is described. 
Systematics affecting previous measurements are addressed by careful characterization of the emis- 
sion spectrum of the neutron source, use of a small scintillator coupled to an ultra-bialkali high 
quantum efficiency photomultiplier, and evaluation of non-linearities in the electron recoil response 
via Compton scattering. A trend towards a rapidly diminishing quenching factor with decreasing 
sodium recoil energy is revealed. Additionally, no evidence for crystal lattice channeling of low- 
energy recoiling ions is found in a scintillator of known crystallographic orientation. A discussion on 
how these findings affect dark matter searches employing NaI[Tl] (e.g., DAMA/LIBRA) is offered. 
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I. INTRODUCTION 

The interpretation of dark matter search data depends 
critically on a complete understanding of detector re- 
sponse to nuclear recoils. These are expected to be in- 
duced via elastic scattering of dark matter particles off 
target nuclei [l[. The modest recoil energies anticipated 
from Weakly Interacting Massive Particle (WIMP) in- 
teractions complicates this necessary labor of detector 
characterization. In recent times, a number of poten- 
tial signals have been reported by several such searches. 
These anomalies [2-4], together with claims of exclusion 
of the relevant WIMP phase space by other experiments 
[il, [q] , have underlined the importance of proper detec- 
tor characterization prior to the reporting of positive or 
negative results. 

Nuclear recoils (NR) induce a diminished response in 
target materials, when compared to electron recoils (ER) 
of the same energy. In scintillators, this is mainly due 
to the small fraction of the recoil energy {Er) that is 
transferred to electron excitation: an energy-dependent 
quenching factor (Q) can be defined as the ratio of the 
light yield from a nuclear recoil to that from an elec- 
tron recoil of the same energy. Nuclear recoil energies 
(typically expressed in units of keVnr) therefore corre- 
late to a smaller observable electron equivalent energy 
(keVee) through this quenching factor. The energy scale 
is typically calibrated using monochromatic gamma rays, 
leading to a satisfactory knowledge of the electron recoil 
response. However, in order to predict the response to 
WIMP-induced nuclear recoils, a dedicated characteriza- 
tion of the quenching factor is necessary. This is typically 
accomplished using neutron-induced nuclear recoils. 

Quenching factor measurements can be divided into 
two groups. The first involves use of sources generating 
a broad spectrum of neutron energies, such as AmBe or 
^^^Cf. A comparison between the simulated and observed 



response to the source is used to generate a best-fitting 
quenching factor, often approximated as being energy- 
independent [Z, iSj- If an attempt is made to include 
the dependence on energy, uncertainties in the simulation 
and systematic effects in the measurement are forced into 
the derived Q{Er). This is specially grav e at the smallest 
energies investigated (few keVnr) [9|, [lQ| . A second more 
reliable method, yet not entirely devoid of its own issues 
[9], involves monochromatic neutron sources. The detec- 
tion of the scattered neutron at a fixed angle from the 
incoming neutron beam leads to a precise knowledge of 
the recoil energy deposited in the target material under 
test [iMi. 
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In the particular case of NaI[Tl] scintillators like those 
used by the DAMA/LIBRA collaboration [2], the first 
method has been reported to yield a sodium recoil 
quenching factor QNa = 0.30 ± 0.01 averaged over 6.5 to 
97 keVnr [7]. For iodine recoils this was Qi = 0.09 ±0.01 
over the range 22 to 330 keVnr [7|. In |8i] the values 
found were QNa = 0.4 ± 0.2 and Qi = 0.05 ± 0.02, av- 
eraged over 5-100 keVnr and 40-300 keVnr, respectively. 
The second method yields values in principle compatible 
with these. However, in the most recent measurement of 
this second type [ll| , a trend towards an increasing QNa 
with decreasing recoil energy is observed down to the 
lowest measured energy, Er =10 keVnr. This increase is 
also predicted by a semi-empirical treatment of the scin- 
tillation mechanisms at play [17]. Indeed, the kinematic 
threshold [18] below which this QNa would be expected 
to decrease should appear at a considerably smaller ~2.5 
keVnr for sodium recoils in NaI[Tl] [19]. In contrast to 
this possibility, other scintillators such as liquid xenon 
(LXe) are expected to display a decreasing quenching 
factor starting already below kinematic thresholds at few 
tens of keVnr [Tij. 

The channeling of recoiling ions within a crystal lattice 
constitutes an additional source of uncertainty specific to 
single crystal scintillators such as NaI[Tl]. If this process 
is present at the few keVnr energies of interest, it can 
result into a considerable displacement of the region of 



WIMP parameter space (mass, coupling) able to explain 
DAM A/LIBRA observations, away from values already 
severely constrained by other experiments |20]. A recent 
theoretical reanalysis of this possibility [21] concluded 
that while ion channeling is a process well-established at 
higher energies, it should play a negligible role for dark 
matter searches. 

The new measurement of QNa and Qi in NaI[Tl] pre- 
sented here attempts to address the uncertainties de- 
scribed above with the introduction of several improve- 
ments, delineated in the next section. These include 
the use of a scintillator with a (partly) known crystal- 
lographic orientation. This allowed a first experimen- 
tal test of the possibility of channeling at few keVnr. 
An increasing Qnr towards low Er and/or non- negligible 
ion channeling would result in a common region of inter- 
est in WIMP phase space for the DAMA, CoGeNT and 
CRESST anomalies [22|, [23|. As discussed below, nei- 
ther of these effects was observed. Based on present re- 
sults, the energy- independent values Qne = 0.3, Qi = 0.1 
that are typically adopted in the interpretation of the 
DAM A/LIBRA annual modulation should be held sus- 
pect. The last section in this paper briefly discusses 
the impact of these findings on the interpretation of 
DAMA/LIBRA results, and their compatibility with 
other searches. 



II. PRESENT IMPROVEMENTS TO THE 
METHOD 

The experimental setup delineated in this and the fol- 
lowing sections was designed to address several sources of 
systematic error affecting previous measurements. These 
measures can be listed: 

o New photomultipliers with enhanced quantum ef- 
ficiency (Hamamatsu's Ultra-Bialkali, UBA) be- 
came commercially available since the last previ- 
ous quenching factor measurement in NaI[Tl] [ll|. 
Their use here results into a light yield of up to 
20 photoelectrons (PE) per keVee, a factor of ~4 
larger than in [lH, i.e., an energy threshold lower 
by the same factor. The reader is referred to the 
discussion around Fig. 4 in [9] for a detailed ex- 
planation on how an insufficient light yield can re- 
sult into artificially large values of Q at the low- 
est recoil energies measured. Taking as a reference 
the Er ~10 keVnr at which the present measure- 
ment starts to suffer from such threshold effects, 
and scaling by the presently achieved increase in 
light yield, it is possible to conclude that all pre- 
vious measurements of Qnu below Er ~40 keVnr 
have been disturbed to some extent by the artifact 
described in [9]. In contrast to this, only the lowest 
Er =8 keVnr data point in this work is expected to 
be somewhat affected. 

o The small dimensions of the NaI[Tl] crystal 
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FIG. 1. Measurement of En using a C-S spectrometer (see 
text). Pulse-shape analysis 31] was unnecessary, in view of 
the clear separation between the ^He(n,p)^H peak and ^He 
recoils. The conventional response to epithermal neutrons 
from a moderated Am/Be source is shown for reference. The 
NaI[Tl] scintillator was positioned in front of the MP320 end 
nearest to the target plane, therefore sustaining a En =2.25 
MeV irradiation. Inset: measurements of neutron yield using 
the rate under the fast-neutron peak in a 1 c.c. ^LiI[Eu] scin- 
tillator [sj. A simulation of the ^Li(n,a)^He reaction rate us- 
ing MCNP converts this measured rate into the isotropic point 
source-equivalent yield depicted, in good agreement with D-D 
generator specifications. Multiple measurements were made 
to monitor the absence of deuterated target depletion during 
neutron scattering runs. Error bars are statistical only, with 
the dominating systematic being the reproducibility in the 
positioning of the ^LiI[Eu] crystal. 



presently used (17x17x27 mm) were sanctioned 
via MCNP-PoliMi Monte Carlo simulation [24] to 
ensure that multiple-scattering played a negligible 
role, comprising just ~10.5% of total neutron in- 
teractions. This is to be compared with the con- 
siderably larger crystals used in previous measure- 
ments (e.g., 5 cm diameter by 5.4 cm long in prj). 
Multiple scattering in oversized detectors creates 
a background continuum able to encumber single- 
scattering signals. This can result in systematic 
effects leading to artificially large values of Q. A 
recent example of how a reduction in detector size 
can help alleviate this background and eliminate 
such systematics can be found in a comparison of 
the results from ^| with those in (26| . 

o The response to low-energy electron recoils in 
NaI[Tl] is known to exhibit a measurable non- 
linearity [27, 28], reaching a maximum in light yield 
at around 10 keVee. While this has been acknowl- 
edged in a few previous quenching factor measure- 
ments [ll|, [13] , this response is typically simply nor- 
malized to that at ~100 keVee, and approximated 
as linear. This action, on its own, results into over- 
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estimating Qnu in the region of energy next to 
DAMA/LIBRA's threshold (-2 keVee) by approxi- 
mately 15% (the precise electron recoil energy scale 
in DAM A/LIBRA can be stablished through a con- 
venient line at 3.2 keVee due to ^^K contamina- 
tion). In this work the response to electron recoils 
is measured throughout the range 2-50 keVee via 
Compton scattering. Emphasis is placed on min- 
imizing alterations of the experimental setup and 
analysis protocol between electron scattering and 
neutron scattering measurements. The quenching 
factor measured here folds in the effect of any non- 
linearities by taking the ratio between light yields 
for NR and ER of the same energy, at all energies 
studied. Unfortunately, this optimal approach to 
extracting Q is not widely used: the extent to which 
a large non-linear response to low-energy electron 
recoils in LXe [29] might affect the interpretation 
of dark matter experiments using that target is an 
open question. 

o Most monochromatic neutron sources exhibit some 
energy dependence on the emission angle, measured 
from the direction of the projectile incident on the 
target sustaining the reaction. In the case of a D-D 
"neutron gun" (deuterium onto a deuterated tar- 



get) like the Thermo MP320 generator employed 
here, this is known to be En ~2.8 MeV in the for- 
ward direction. En -2.4 MeV at 90°, and En --2.2 
MeV at 180°, for an accelerating potential of 80 kV 
[30] . Information about the sense of the deuterium 
projectile along the longitudinal axis of the gun was 
however not available from the manufacturer. The 
neutron energy emitted towards the NaI[Tl] scin- 
tillator ("axial" in Fig. 1) was measured using a 
large- volume ^He counter (9.2 cm diameter, 31 cm 
long). A counter this large is able to contain the 
full projected range of protons generated by the 
^He(n,p)^II reaction. This allows for low-efficiency, 
high-resolution spectroscopy of monochromatic fast 
neutron sources, for which a full-energy deposition 
peak is detected at En + 764 keV, where 764 keV 
is the Q-value of this exothermic reaction [3l| (Fig. 
1). ^He counters used for this unconventional ap- 
plication are referred to as Cuttler-Shalev (C-S) 
spectrometers [31]. This characterization allowed 
to remove an uncertainty in En that may have af- 
fected previous use of the Thermo MP320 [32, 
Additional measurements of the stability of neu- 
tron yield were performed using a small ^LiI[Eu] 
scintillator (Fig. 1). 
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FIG. 2. Experimental setup used for Compton scattering measurements (see text). 



III. COMPTON SCATTERING 
MEASUREMENT 



Fig. 2 displays the setup used to measure the re- 
sponse of NaI[Tl] to low-energy electron recoils. A pen- 
cil beam (0.6° aperture) of collimated 356 keV gam- 
mas from a 0.1 mCi ^^^Ba source bathed the NaI[Tl] 
scintillator. Gamma source and lead collimator were 
mounted onto a goniometric stage used to select the 
Compton scattering angle, Be- The square cross-section 



(17x17x27 mm) Czochralski-grown NaI[Tl] single crys- 
tal was obtained from Proteus/ Amcrys. It contained ap- 
proximately 700 ppm of thallium dopant. This is simi- 
lar to DAMA/LIBRA's quoted -0.1% [35] (smah vari- 
ations in dopant concentration are not expected to af- 
fect quenching factor measurements [36]). A request was 
made to the manufacturer to have two of its side faces be 
cut perpendicular to the [100] crystal growth axis (Fig. 
5). The scintillator was vertically mounted onto a Hama- 
matsu ultra-bialkali (UBA) R7600U-200 photomultiplier 
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(PMT), using optical RTV as a couplant. PMT and base 
were mounted onto a second goniometric stage ^] used 
to rotate them around their vertical during ion channel- 
ing measurements (Fig. 5, section VI). Both were sur- 
rounded by a sheath of mu-metal to avoid magnetic field 
disturbances to the PMT gain upon rotation. The preci- 
sion of the angular alignment of all elements involved is 
estimated at ~ 1°. 

A second lead collimator (1.7° aperture) was placed in 
front of a 53.3 mm diameter, 40.5 mm long intrinsic ger- 
manium (HPGe) detector, used to capture the Compton- 
scattered gamma. Careful alignment of both collimators 
was performed at 6c = 0°. Fig. 3 contrasts the energy 
loss to electrons in NaI[Tl] with its Klein-Nishina expec- 
tation, as a function of 9c. 
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FIG. 3. Difference between ^^^Ba incident gamma energy 
(356 keV) and the centroid of the distribution of HPGe sig- 
nals coincident with NaI[Tl] events (i.e., energy deposited in 
NaI[Tl] via Compton scattering), plotted against Be- Er- 
ror bars are not statistical: horizontal bars correspond to the 
aperture of the second lead collimator, and vertical to the 
standard deviation in the HPGe signals. Their overlap with 
the Klein-Nishina expectation is witness to the good angular 
alignment obtained in this setup. Inset: Example spectra of 
HPGe signals coincident with NaI[Tl]. Two peaks used for 
HPGe energy calibrations are shown. These were performed 
prior to each Gc measurement, to monitor HPGe gain stabil- 
ity. 



action, and compared to the single photoelectron (PE) 
mean current (Fig. 8 inset), in order to express the 
NaI[Tl] energy deposition in number of PE. The electron 
recoil energy deposited in NaI[Tl] was obtained from the 
difference between the incident 356 keV gamma energy 
and the energy registered in the HPGe (horizontal axis 
in Fig. 4). The inset in Fig. 4 displays the normaliza- 
tion of these results to PE yield per keVee imparted to 
an electron recoil. The expected nonlinearity [28| 
is observed. The quantum efficiency of the UBA PMT 
is estimated at ~33% for the NaI[Tl] spectrum of light 
emission: accounting for the ~25% reduction in light 
yield expected for 1 MeV gammas this result is in 
rough agreement with expectations (47 scintillation pho- 
tons/keVee here vs. 40 photons/keVee typically quoted 
for NaI[Tl] at 1 MeV). Given the good agreement of the 
observed PE yield with the corresponding curve in 
(Fig. 4, inset), this last is adopted to extrapolate present 
measurements to the interval 50 - 200 keVee, via the ex- 
pression PE/keVee = 16.65 + 3.81 x e-O-Oi^sx/ceVee^ 
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FIG. 4. Number of PE generated by Compton-scattered elec- 
trons in NaI[Tl] vs. energy registered in HPGe (see text). 
Inset: NaI[Tl] photoelectrons per keVee measured in the 2-50 
keVee range. Red dots correspond to the mean value com- 
puted for 2 keVee bins. The red line is the non-linear behavior 
described in [27], normalized to these data at 10 keVee. 



Binary data from the Acqiris digitizer were read and 
analyzed using a Lab VIEW code. Fig. 4 shows the sum 
of all Compton measurements. A band of coincident 
IIPGe-NaI[Tl] signals is clearly discernible against resid- 
ual backgrounds, following the application of cuts to the 
data (removal of events with significant NaI[Tl] scintil- 
lation present already prior to coincidence time, and of 
NaI[Tl] signals with mean scintillation decay times out- 
side of the 50-500 ns range). The current at the NaI[Tl] 
PMT was integrated over the 3 jas following an inter- 



IV. NEUTRON SCATTERING MEASUREMENT 

The arrangement in Fig. 2 is converted to that in Fig. 5 
by simple removal of HPGe detector, gamma source and 
collimators, and switching of two cables. The same gonio- 
metric stage used to hold the gamma source is utilized to 
mount a Bicron 501A liquid scintillator cell (5 cm diame- 
ter, 5 cm long), able to distinguish neutron from gamma 
interactions via pulse-shape discrimination (PSD) [38] 
(Fig. 6). This cell is used to detect neutrons scattered 
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off NaI[Tl] at an angle It is also optimal in its decay 
time characteristics (3.2 ns) for fast-coincidence applica- 
tions. Its PMT was surrounded by mu-metal: tests were 
performed to ensure that no measurable change in PMT 
gain occurred with A large tank of borated water 
was used to block D-D neutrons from reaching the cell 



following a direct path. The distance between 501A cell 
and NaI[Tl] was varied within the interval 20 - 45 cm, 
using the shortest distances for large values of <l>n, in or- 
der to compensate for the smaller number of scattered 
neutrons expected there, at the expense of some loss in 
angular resolution. 
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FIG. 5. Experimental setup used for neutron scattering measurements (see text). The long (vertical) axis of the NaI[Tl] crystal 
points into the plane of the figure. A thin lead foil, not shown, is used to block low-energy x-rays emanated by the D-D 
generator. The distance between deuterated target plane and NaI[Tl] was 42 cm. The PS776 is a 16 channel amplifier. 



An important detail noticeable in Fig. 5 is the trigger 
configuration. The coincidence sequence between 501A 
and NaI[Tl] detectors is armed by events in the first, and 
closed by events in the second, i.e., in time-reversed or- 
der. To compensate for this, the NaI[Tl] signal is delayed 
by ~250 ns, and the TAC/SCA used to determine coinci- 
dences is set to a sufficiently long window (500 ns). Digi- 
tizer trace length and external trigger position within are 
set to capture the entire sequence of events. The reason 
for this counterintuitive arrangement is the long decay 
time of NaI[Tl] scintillation (~200 ns), when compared 
to the time-of-flight (TOF) for -2.2 MeV scattered neu- 
trons traveling between the two detectors (—15 ns. Fig. 
7). A standard trigger arrangement where NaI[Tl] events 
arm the TAC/SCA results in many missed triggers for 
low-energy (few PE) NaI[Tl] signals where no PE is gen- 
erated within the first 15 ns following a nuclear recoil. 
This systematic effect is entirely bypassed by the present 
arrangement. 

NaI[Tl] PMT current integration is performed identi- 
cally for Compton scattering data, single PE spectrum 
reconstruction (Fig. 8, inset), and neutron scattering 
runs. This is important, given that the quenching factor 
is determined by the ratio of these currents (expressed 
in number of PE equivalent) for comparable electron and 
neutron recoil energies. To this end, any small fluctu- 
ations in the DC-level of NaI[Tl] PMT traces are ac- 
counted for in the analysis, PMT current integration 
times are kept constant for NR and ER, etc. The ex- 
perimental arrangement emphasizes, by design, keeping 
changes to a minimum in going from electron scattering 
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FIG. 6. PSD against gamma backgrounds affecting the Bicron 
501 A cell 33. 38]. Events below the dotted line are passed 
by the off-line analysis, maximizing the acceptance of scat- 
tered neutrons at the expense of some gamma contamination. 
More stringent cuts (dashed line) are used for measurements 
of CsI[Na] quenching in this same setup [39| . 



to neutron scattering measurements. Background cuts 
during the analysis of neutron recoil data are the same 
as in section III, with the addition of those described in 
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FIG. 7. Events passing all cuts in neutron scattering runs 
taken at = 44°, depositing an average 54 keVnr in NaI[Tl]. 
Inset: distribution of events in the 10 < PE < 400 interval, 
projected onto the horizontal time axis. The expected TOP 
of 15 ns for a En = 2.25 MeV neutron traversing a distance 
d = 30 cm between NaI[Tl] and 501 A detectors is observed 
[ll[. See text for a discussion on trigger configuration. 



Fig. 6. The resulting events passing cuts display a clear 
separation between nuclear recoils and residual NaI[Tl] 
- 501A coincidences mediated by numerous thermal cap- 
ture background gammas (Fig. 7). 



V. EXTRACTION OF THE QUENCHING 
FACTOR 



Three vertical dotted lines in the inset of Fig. 7 define 
two time domains of equivalent span, the one on the left 
containing a background of spurious events, the one on 
the right dominated by nuclear recoil signals in NaI[Tl]. 
A spectrum of NaI[Tl] light yield (number of PE) can be 
formed for events in each of these time domains, for each 
scattering angle probed. Subtraction of the background 
spectrum from the recoil-dominated spectrum generates 
the residuals in Fig. 8. The centroids of Gaussian fits 
to these provide the mean number of PE generated by 
nuclear recoils in NaI[Tl] at each The correspond- 
ing mean nuclear recoil energies (horizontal axis values in 
Fig. 9) are extracted from MCNP-PoliMi [2i] simulations 
of the runs. These energies are in good agreement with 
expectations from basic neutron kinematics [40]. The 
ratio between these PE centroids and the PE yield ob- 
served for electron recoils of the same energy (Fig. 4, and 
extrapolation beyond 50 keVee discussed in section III) is 
identified as the sought quenching factor. This is shown 
in Fig. 9. 

Fig. 8 also displays the calculated triggering efficiency 
as a function of number of PE generated by a neutron 
recoil in NaI[Tl]. Formally, this is computed as the prob- 
ability of a binomial distribution with PE trials and com- 
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pound success probability p = piP2 returning a number 
of successes equal or larger than one, where pi = 0.69 is 
the probability of an individual PE triggering the CFD 
(Fig. 8, inset) and P2 = I - exp{-At/r). Here At = 230 
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ns is the time window available for a PE to set off the 
trigger configuration described above (see Fig. 7) and r ~ 
200 ns is the decay time for NR-induced scintillation (Fig. 
10). In other words, p2 accounts for the finite probability 
that the first PE generated in the NaI[Tl] is not prompt 
enough to stop the TAC/SCA within its 500 ns range, 
after accounting for the effect of the delay loop. The ef- 
fect of p2 on reducing the triggering efficiency for a given 
number of PE is more pronounced for slower (longer r) 
scintillators such as CsI[Na] (39| . 

As computed, this triggering efficiency reaches unity 
for PE > 10, i.e., only the PE spectrum corresponding 
to runs at the smallest scattering angle ^n= 16° {Er'^S 
keVnr) requires a correction based on it. That the rest of 
the measurements collected at larger scattering angles are 
free from any such threshold effects can be ascertained by 
observing the left shoulders of the distributions in Fig. 8, 
which are not similarly shaped by a triggering efficiency 
curve (see discussion around Fig. 4 in [9]). 

Iodine recoils emerged distinctly above the few PE spu- 
rious coincidence noise (Fig. 7) only at the three largest 
values of attempted (72°, 78°, 102°). Given the 
proximity of their PE distributions to the triggering effi- 
ciency threshold, these were all corrected by the efficiency 
curve prior to extracting a Gaussian best-fit. The iodine 
quenching factor obtained is considerably smaller than 
for sodium recoils, similarly to previous results 0, Q, 
and comparable to values obtained for Cs and I recoils 
in CsI[Na] in this same setup [39]. 
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VI. SEARCH FOR AN ION CHANNELING 
EFFECT 

The goniometric stage holding the NaI[Tl] detector 
can be used to select an orientation allowing a recoil- 
ing sodium ion to (in principle) channel down the (100) 
plane perpendicular to the known [100] growth axis of the 
crystal (Fig. 5 and inset of Fig. 11). The corresponding 
angle ^rec varies with <l>^, and can be calculated using 
basic kinematic relations [40]. In [20^ it was postulated 
that such a channeled ion should exhibit a QNa i-e., 
lose energy exclusively through ionization much like an 
electron, and not via secondary nuclear recoils. In [2l| 
it is claimed that "blocking" effects should render this 
impossible for recoiling ions originating in nuclei initially 
at rest on the lattice, i.e., the case for WIMP or neutron 
interactions. 
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FIG. 11. Results from a negative search for sodium recoil 
channeling in single-crystal NaI[Tl] (see text). "Random ori- 
entation" corresponds to an alignment of the [100] crystal axis 
with the direction of the incoming neutrons, i.e., an orienta- 
tion for which no significant channelin g w ould be expected. 
The inset displays a HEALPix [2l|, [4l,T43[| map of available 
channels in the NaI[Tl] lattice, represent by red bands. These 
are to scale for Er — 10 keVnr (see text). 



FIG. 10. Grayscale intensity plots showing the observed mean 
decay time of NaI[Tl] scintillation. NR energies are in keVnr, 
ER in keVee. A single exponential decay component is used 
in fits to the running integral of PMT currents. The bottom 
panel shows the centroids and dispersion of the data in 2 keV 
bins. The evolution of these centroids with energy is very 
similar to that found in previous work [li,ji,[il,gl|. 



The angular precision of the present setup cannot be 
claimed to be better than ~ 1°. Similarly, the manu- 
facturer of the NaI[Tl] scintillator listed the tolerance of 
the requested crystal orientation to be "a few degrees". 
Nevertheless, the angular dimension of the (100) planar 
channel would be sizable for low-energy sodium recoils. 
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For example, in NaI[Tl] this room-temperature channel- 
ing angle should be ~ ±5.6° around the (100) plane for 
Er = 10 keVnr, - ±4.0° at Er = 28 keVnr, and - ±3.5° 
at Er = 4:l keVnr [42]. Taking into consideration the ef- 
fect of the finite solid angle subtended by the 501A cell 
as seen from the position of the NaI[Tl] scintillator, the 
measurement attempted here should be forgiving enough 
to display some evidence for ion channeling, if the pro- 
cess was possible at few keVnr for nuclei knocked off their 
positions in the lattice. 

Fig. 11 illustrates the negative results from this search. 
Vertical arrows mark the expected position of the light 
yield distribution for sodium recoils if channeling was re- 
alized, i.e., for QNa = 1- No difference was observed 
between control runs ("random orientation") and runs 
with values of ^rec selected to favor channeling along 
(100). For = 28° additional runs at ^^ec = 72° and 
^rec = 78°, i.e., a few degrees beyond the nominally fa- 
vored ^rec = 75°, were performed. This was done to 
account for the possibility of some significant misalign- 
ment in the setup. These runs returned spectra similar 
to those shown in the figure. 

Following these measurements, an attempt was made 
to test the crystallographic alignment of the crystal using 
x-ray scattering at Argonne National Laboratory. The 
presence of Laue spots confirmed the single-crystal na- 
ture of the sample, but it was impossible to obtain in- 
formation beyond the few degree tolerance listed by the 
manufacturer (x-ray penetration through the scintillator 
casing was minimal, and its removal led to the expected 
rapid degradation of hygroscopic NaT [Tl] surfaces). How- 
ever, for the reasons listed above, and barring a neglect 
by the crystal manufacturer to adhere to alignment in- 
structions, the absence of any excess at QNa = ^ in the 
red histograms of Fig. 11 suggests that the arguments 
against channeling in [2l| are now experimentally con- 
firmed. 



VII. IMPLICATIONS FOR DARK MATTER 
SEARCHES 

While the decreasing Q Na obtained in this experiment 
challenges several previous results, the author has con- 
firmed a value of Qatc ^ 0-1 for recoils below 24 keVnr 
using a new independent technique involving an ^^Y/Be 
photo-neutron source. This is treated in a separate publi- 
cation [i^ . As discussed in section II and [9] , the combi- 
nation of poor light yield and lack of sufficient control of 
systematics near threshold can produce a false impression 
of constant or increasing quenching factors with decreas- 
ing recoil energy. Based on present experimentation, this 
seems to have been the case for previous measurements 
of Qnq below Er ~ 40 keVnr, from which a discussion of 
threshold effects is notoriously absent (specifically, con- 
trol of expected vs. observed event rates at recoil energies 
strongly impacted by triggering efficiency [9]). In recent 
times, improvements to the methodology of quenching 



factor measurements for LXe [25|, |26| have clarified a sim- 
ilar situation, one that led to an interpretation of dark 
matter search data [45( based on markedly optimistic 
quenching factors [195. New quenching measurements 
performed by the author on CsI[Na] [39] indicate that 
this unfortunate situation may also extend to previous 
studies of this other scintillator i36ll . 
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FIG. 12. Effect of QNa on a DAMA/LIBRA ROI, adapted 
from ^2S\. XENONIOO exclusions have present limited credi- 
bility for < 12 GeV/c^ A CDMS ROI is from [5l|. 

Fig. 12 encapsulates the impact of present measure- 
ments specifically on DAMA/LIBRA. A value of Qj^a ^ 
0.1 for putative WIMP recoils leading to energy depo- 
sitions near the 2 keVee DAMA/LIBRA threshold dis- 
places the region of interest (ROI) for a WIMP in- 
terpretation well into the realm excluded by several 
other techniques, significantly increasing the existing 
tension. While deviations from the assumption of a 
Maxwellian velocity distribution for the local galactic 
halo can still reconcile this displaced DAMA/LIBRA ROI 
with CRESST and CoGeNT anomalies @, 0, this is 
made more difficult by present results. It can be argued 
that these halo velocity deviations are to be expected 
[48] : the upcoming availability of an improved catalog of 
nearby stellar velocities from the GAIA satellite probe 
should be able to remove these uncertainties (49| . 

The effect of ion channeling would have been to dis- 
place the DAMA/LIBRA ROI towards smaller values of 
WIMP mass and lower couplings [20|, i.e., precisely to- 
wards agreement with other anomalies, and away from 
constraints from other techniques. The objections raised 
in [21] and the present absence of any indication for this 
process leave little room for this possibility (notice how- 
ever the possible caveat listed above). 

In view of the state of affairs, our emphasis should 
be on continuing to improve the quality of these criti- 
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cal low-energy calibrations [44^. Efforts [46*] to interpret 
the microscopic physics underlaying existing (but ques- 
tionable) calibrations may be premature. If the trend 
towards vanishing quenching factors at few keVnr contin- 
ues, scintillating materials may prove to be of limited to 
no use in searches for low-mass (^10 GeV/c^) WIMPs. 

ACKNOWLEDGMENTS 

The author is indebted to N. Fields, D. Hooper, T. 
Hossbach and N. Weiner for many useful conversations. 



N. Bozorgnia, G. Gelmini and P. Gondolo for encourage- 
ment and much important information about channeling, 
and F. Bartolome, M.L. Sarsa and X. Huang for sug- 
gesting and facilitating x-ray scattering measurements 
at ANL. The author acknowledges the Aspen Center for 
Physics, supported by NSF Grant 1066293, for hospital- 
ity during the completion of this manuscript. 



[1] M.W. Goodman and E. Witten, Phys. Rev. D31 (1985) 
3059; I. Wasserman, Phys. Rev. D33 (1986) 2071. 

[2] R. Bernabei et a/., Eur. Phys. J. C56 (2008) 333. 

[3] C.E. Aalseth et a/., Phys. Rev. Lett. 101 (2008) 251301; 
Erratum ibid 102 (2009) 109903; Phys. Rev. Lett. 106 

(2011) 131301; Phys. Rev. Lett. 107 (2011) 141301. 

[4] G. Angloher et al, Eur. Phys. J. C72 (2012) 1971. 

[5] E. Aprile et a/., Phys. Rev. Lett. 109 (2012) 181301. 

[6] Z. Ahmed et a/., Phys. Rev. Lett. 106 (2011) 131302. 

[7] R. Bernabei et al, Phys. Lett. B389 (1996) 757; R. Bern- 
abei et a/., Riv. Nuovo Cimento 26 (2003) 1. 

[8] K. Fushimi et al, Phys. Rev. C47 (1993) 425. 

[9] J.L Collar, arXiv: 1010.5187, 
[10] J.L Collar. arXiv : 1106.0653. 
[11] H. Chagani et al, JINST 3, P06003 (2008). 
[12] Th. Jagemann et a/., Nucl. Instr. Meth. A564 (2006) 
549. 

[13] E. Simon et a/., Nucl. Instr. Meth. A507 (2003) 643. 

[14] C. Gerbier et al, Astropart. Phys. 11 (1999) 287. 

[15] D.R. Tovey et a/., Phys. Lett. B433 (1998) 150. 

[16] N.J.C. Spooner et a/., Phys. Lett. B321 (1994) 156. 

[17] V. L Tretyak, Astropart. Phys. 33 (2010) 40. 

[18] D.J. Ficenec et a/., Phys. Rev. D36 (1987) R311; S.P 

Ahlen et al, Phys. Rev. D27 (1983) 688. 
[19] J.L Collar and D.N. McKinsey, arXiv: 1005 . 0838 
[20] R. Bernabei et a/., Eur. Phys. J. C53 (2008) 205; E.M. 

Drobyshevski, Mod. Phys. Lett. A23 (2008) 3077. 
[21] N. Bozorgnia, G.B. Gelmini and P. Gondolo, JCAP 11 

(2010) 019. 

[22] D. Hooper et a/., Phys. Rev. D82 (2010) 123509. 

[23] C. Kelso, D. Hooper, M.R. Buckley Phys. Rev. D85 

(2012) 043515. 

[24] S.A. Pozzi et al, Nucl. Instr. Meth. A513 (2003) 550. 
[25] E. Aprile et al, Phys. Rev. C79 (2009) 045807; A.G. 

Manalaysay, Ph.D. Diss., University of Florida, 2009. 
[26] G. Plante et a/., Phys. Rev. C84 (2011) 045805. 
[27] W.W. Moses, Nucl. Instr. Meth. A487 (2002) 123. 
[28] I.V. Khodyuk et a/., J. Appl. Phys. 107 (2010) 113513; 

W.S. Choong et a/., IEEE TNS 55 (2008) 1073; W.W. 



Moses et a/., IEEE TNS 55 (2008) 1049; L.R. Wayne et 
a/., Nucl. Instr. Meth. A411 (1998) 351; K. Saito et a/., 
Nucl. Instr. Meth. 185 (1981) 299; F.T. Porter et a/., 
Nucl. Instr. Meth. 39 (1966) 35. 
[29] E. Aprile et al, arXiv: 1209.3658 

[30] "CRC Handbook of Fast Neutron Generators", Vol. 1, 
by J. Csikai. CRC Press, Boca Raton, FL, 1987. 

[31] A. Sayres and M. Coppola, Rev. Sci. Instr. 35 (1964) 431; 
D.L. Chichester, J.T. Johnson and E.H. Seabury, 21st 
Intl. Conf. on the Application of Accelerators in Research 
and Industry Fort Worth, TX, August 2010. 

[32] J.A. Nikkei et al, Astropart. Phys. 29 (2008) 161. 

[33] A. Manzur et al, Phys. Rev. C81 (2010) 025808. 

[34] R.B. Murray Nucl. Instr. 2 (1958) 237. 

[35] R. Bernabei et al, Nucl. Instr. Meth. A592 (2008) 297. 

[36] H. Park et al, Nucl. Instr. Meth. A491 (2002) 460. 

[37] Model RSP-IT, Newport Corp, Irvine, CA. 

[38] D. Wolski et al, Nucl. Instr. Meth. A360 (1995) 584. 

[39] CosI collaboration, in preparation. 

[40] Y. Shimizu et al, Nucl. Instr. Meth. A496 (2003) 347. 

[41] PF. Smith et al, Phys. Lett. B379 (1996) 299. 

[42] N. Bozorgnia, priv. comm. 

[43] http : //healpix . j pi . nasa . gov/ 

[44] J.L Collar, arXiv: 1303.2686, submitted to Phys. Rev. 
Lett. 

[45] E. Aprile et al, Phys. Rev. Lett. 105 (2010) 131302. 

[46] P. Sorensen, JCAP 10 (2010) 033; P. Sorensen and E. 
Dahl, Phys. Rev. D83 (2011) 063501; M. Szydagis et al, 
JINST 6 (2011) 10002. 

[47] P.J. Fox et al, Phys. Rev. D85 (2012) 036008; P. Gon- 
dolo and G. Gelmini, JCAP 12 (2012) 015. 

[48] M. Kuhlen, M. Lisanti and D.N. Spergel, Phys. Rev. D86 
(2012) 063505; C.W. Purcell, A.R. Zentner and M.Y. 
Wang, JCAP 08 (2012) 027; K. Freese, M. Lisanti and 
C. Savage, Rev. Mod. Phys., in press, arXiv: 1209.3339 

[49] C. Frenk, priv. comm. 

[50] M. Dress and G. Gerbier, PDG 2012 dark matter review, 

arXiv: 1204.2373 
[51] J.L Collar and N.E. Fields, arXiv: 1204.355. 



